Dark states of dressed spinor Bose-Einstein condensates of spin-1 atoms 
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Dressed states of spinor Bose-Einstein condensates of spin-1 atoms located in optical cavity of 
near resonance frequency of excited state of the atom are investigated. The wave function of dark 
states are given. 
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Implementation of optical trap of spin-1 to spin-3 
bosonic atoms, which librate the freedom of spin of the 
atoms, has given a great chance to investigate more 
prolific physics of cold atomic cluster p], 0, H, 0, @. 
When temperature is low, spinor Bose-Einstein conden- 
sate (spinor BEC) is form, whose eigenstate conserve the 
total spin and its Z-component [6| . If the system is cou- 
pled to a high Q-factor optical cavity mode, the eigen- 
state will be modified. This is a system of cavity quantum 
electrodynamics (cavity QED) with a spinor BEC. 

Cavity QED investigate the interaction between atom 
and optical cavity mode 0, H[ . In the parameter region 
that the interaction strength between atom and optical 
mode is much larger than decay rate of cavity mode and 
excited state of the atom, the system is in strong coupling 
region. In this region, energy transfer between atom and 
cavity mode by absorbing and emitting photon periodi- 
cally. The eigen state of the coupling system form dressed 
state described by Jaynes-Cummings model |9|. When 
there are more than one atom, the coupling effect will 
be enhanced. Several previous research has investigated 
cavit y Q ED with a BEC theoretically and experimen- 
tally dfl, EMI El, El El El] • Dressed state of BEC is 
confirmed |l5j. For a multicomponent condensates, it is 
shown that dark states of dressed BEC exist [16[ . When 
the dressed BEC is in dark state, no atom is state at 
excited state, but oscillate between two ground states of 
A type energy level atom. The research of dressed BEC 
with has potential in understanding of quantum state fl7j 
and quantum information. 

This paper consider the dark states of dressed spinor 
BEC of spin-1 atoms. At first, coupling Hamiltonian be- 
tween spinor BEC and optical cavity mode is given. And 
then, we calculate the matrix elements of the coupling 
Hamiltonian, where we make use of fractional parentage 
coefficients (FPCs) to get a simple analytic expression. 
At last, wave function of the dark states is deduced and 
value scope of quantum number for the dark states is fix. 

The system we are going to investigate consist of N 
87 Kb atoms optically trapped to form a spinor BEC, lo- 
cated inside a high-Q optical cavity. In the strong cou- 



pling region, the photon decay rate of the cavity and 
the non-resonant decay rate of the atoms is neglected. 
The optical trap can be modeled by a harmonic trap 
potential. The frequency of the cavity mode is near res- 
onance to the transition frequency of D% line of &7 Rb : 
5 2 Si/ 2 ,F = 1 — ► 5 2 P 3/2 ,F = 0, so that only the two 
levels is considered in the dressed states. The cavity con- 
sist of two parallel high reflection mirror, which has two 
optical modes at the resonance frequency: one mode is 
left circular polarization, the other mode is right circular 
polarization. The two modes have opposite angular mo- 
mentum along the axis of the cavity. The Hamiltonian of 
the atoms and cavity system can be written, respective, 
as: 
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where Pi is momentum of the ith atom, m is mass of 
the atoms, U(r) is harmonic potential of optical trap, 
|ej) is the excited state 5 2 P 3 / 2 ,F — of the ith atom, 
huj e is energy level of the excited state, V(ri — r/) is in- 
teraction between two atom that depend on coupling of 
their spins, and hw c is energy level of cavity mode, 
is annihilate operator of standing wave cavity mode of 
polarization fi = ±. V(ri — Tj) can be given as V gg (ri — 



r j)\9i9j){9i9j\ + V ge (ri - r j ){\g i e j ){g i e j 
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where \gi) is ground state of the ith atom. V gg is spin de- 
pendent interaction between two ground state atoms with 
spin-1, which is (en + C2F,; • Fj)<5(rj — r^), Fj is spin oper- 
ator of the ith atom, 00(02) is spin independent (depen- 
dent) interaction constant [l|. For V ge , one of the atom 
is spin-0, the scattering has only one channel, so that the 
interaction is spin independent, i.e. V ge = C4<5(r^ — r^) 
where C4 is a constant. 

The direction of quantize axis of the atoms is along the 
axis of the cavity, so that by absorbing a left (right) circu- 
lar polarization photon from the cavity, only the a + ((t~) 
transition that conserve the angular momentum could 
occur; the same conservation law of angular momentum 
apply for emission one photon to the cavity. Therefore, 
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the Hamiltonian of coupling between atoms and cavity 
mode is, 

N 

H AC = ^2[m + \e i )(g i> -l\a + +m-\e i )(g h +l\a^]+H.c. 

i=l 

(3) 

where Ml± — d±y / huj c /2eoV is coupling strength of 
transition, d± is the corresponding atomic dipole-matrix 
element that satisfy d + = — d_, V is mode volume of the 
cavity mode, |(^,±i) is the ground state of atom with 
spin component ±1. 



The eigenstate of non-coupling Hamiltonian 
Hj\ + Hq can be labeled as \g, S, M; N 1} N 2 ) and 
|e, S', M'\ Ni, N%), where g means all N atoms are at 
ground state, e means one of the atoms is at excited 
state, S and M is total spin and the Z component, Ni 
{N 2 ) is photon number of /i = + (/* = —) polarization 
photon. For the state that all N atoms are at ground 
state, N — S is limited to be even, and the eigen energy is 
approximated to be E^ 2 = J g S(S + l) + ?w c {Ni + N 2 ) 
@, where J g is spinor interaction constant that depend 
on total spin and spatial wave function. When one atom 
is at excited state of spin-0, all the other atoms couple 
to total spin S', with the limitation that N — 1 — S' is 
even. The atomic part of wave function of this state can 
be written as 

\e,S',M') = ^£{b e (rOn^)]<;M'®} W 

where f e (g) is spatial wave function of excited (ground) 
state atom, $g~Mi (?) is total symmetric spin state of all 
the other atoms without the ith atom with quantum num- 
ber 5" and M'. Spin of the ith atom is zero, and it cou- 
ple with spin state of the other N — 1 atoms, forming 
total spin state S'. Using variation method to deduce 
Schrodinger equation for the spatial wave function ip e i g ) 
and solving the equation by perturbation approximation, 
eigen energy of |e, S', M') plus energy of the cavity mode 
is given as E^ 2 = J e S'(S' + 1) + hw e + fi^ c (iVi + N 2 ). 
J e is spinor interaction constant. When N is small, J g 
and J e is approximated to be equated to each other, 
noted as J. J is much smaller than ftw c and ftiu e , so 
that for a fix set of quantum number (<?(e), JVi, N2), the 
energy levels expand to be a narrow band by quantum 
number S(S') around energy level (ftus e ) + huj c (Ni + N 2 ). 
When there is no static magnetic field, the energy levels 
do not depend on M(M'). We consider the dressed states 
consist of states in the first excited band: \g, S, M; 1, 0), 
\g,S,M;0,l) and \e,S',M'\ 0,0). According to the con- 
servation of angular momentum at quantize axis, the 
dressed states with total angular momentum along quan- 



tize axis Mh can be written as 

\M) = ^ S e'' M (t)\^S',M;0,0) 

S' 

+ J2c S + M - 1 (t)\g,S,M-l;l,0) 
s 

+ J2c s J M+1 {t)\g,S,M + 1;0,1) (5) 
s 

where the summation of cf ,M run through S' that 
N — S' — 1 is even and S' >= M, the summation of 
c^' M 1 (c 5 ^' M+1 ) run through S that N — S is even and 
£?>=M-l(S>=Af + l). 

In order to get the equation of motion of probabil- 
ity amplitudes under the full Hamiltonian H = Ha + 
Hc + Hacj we deduce the matrix element of the coupling 
part of Hamiltonian Hac a t the bases of eigen states of 
non-coupling Hamiltonian. We consider matrix element 
(e, S', M'; 0, 0\H AC \g, S, M ; 1, 0) first. The atomic state 
(e,S,M\ is expanded as equation (j4|). The atomic spin 
state of \g, S, M; 1,0) is total symmetric spin state of N 
atoms, i?g M . Making use of FPCs, the spin state of the 
ith atom can be extracted, so that the atomic state can 
be expanded as 

JV 

| 5 , S, M) = [[] <p g ( rj )] [o£ [ X (i)^l M (})}s 
3=1 

+ b%\xWs-iM®]s] (6) 

where af = {[1 + (-l) w - s ](iV - S)(S + 1)/[2JV(2S + 
l)]} 1 / 2 and b% = {[l + {-l) N - s ]S(N + S+l)/[2N(2S + 
l)]} 1 / 2 is FPCs given in reference [3, x(i) is spin state of 
the ith spin-1 ground state atom. The coupling between 
X(i) and $£±1 M (i) can be expanded by using Clebsch- 
Gordan coefficients [l9| . When we calculate the matrix 
element of operator of the ith atom in H A Ci we can ex- 
tracted spin state of the ith atom and make use of the 
orthogonal relation of spin states to get the exact expres- 
sion of the spin part of the element. Under approxima- 
tion that the spatial wave functions are independent of 
total spin and internal state, the matrix element can be 
expressed as 

(e,S',M';0,0\H AC \g,S,M;l,Q) 

+ VNm+bgC§?f liM ,,_ 1 5s,s>+i8M,M>-i (7) 

where C§ M g 2 M Ma is Clebsch-Gordan coefficients. Using 
similar analysis, the matrix element with another photon 
state can be obtained as 

(e, S', M'\ 0, 0\H AC \g, S, M; 0, 1) 
= V r /V?i0_a5 C§Jj_ jM ,a5s,s>-i8m,M'+i 
+ VNm-b%C§% iM , tl 8s,s>+i8M,M>+i (8) 
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Inserting the state into Schrodinger equation, and take 
the inner product with all bases in equation (4), the equa- 
tion of motion of probability amplitudes are obtained as 



ih 



dc 



S',M 



(It 



= [tVjJ e + JS'(S' + 1)] 



JS'.M 



S'-lM-l S'-l,M-l 



S'+1,M—1 S'+1,M—1 



ATkn N nS'-\,M+\ S'-1,M+1 

Ivno h N n s ' +1 ' m 

iV fi-"-»s' + l°S',l,M,l 

where N — S' — 1 is limited to be even, and 



(9) 



S.M-l 



(It 



= [hw c + JS(S +l)]c^ M ~ 1 



v - v u ',, '•" (io) 



ih- 



dc 



S.M+l 



dt 



= [fyjj c + JS(S + 1)] 



S.M+l 



Nm-a%Cg!£l 



+i 



S-l.M 



(11) 



where A — 5 is limited to be even. Mh is total angular 
momentum along quantize axis. The exist of dark state 
require that the probability amplitude cf ' M equate to 
zero. From equation (jTQJ) and (fTTj) , this will result in 
the total spin dependence of the time factor of c£' M=Fl . 
From equation ([9]), this will result in relations between 
different probability amplitude. Because the time factor 
of c|' M=Fl depend on S as e ~^S(S+i)t/h ) and n+ _ _ fi ^ 
equation ([9]) reduce to two equation for dark states, 

N n S'-l,M-l S'-l.M-l AT f-tS'-l.M+l S'-l.M+l _ 



jv — 1 .III — J 

a S'-l°S' 1,M,-1 



t AT ^S'+l.M-l S'+1,M-1 .JV n S' + l,M+l S'+1,M+1 _ n 
°S'+l°S',l,Af,-l c + °S"+1°S',1,M,1 c - u 



(12) 
+i _ 
(13) 



It is obvious that one kind solution is M = and S 1 > 0. 
When S = 0, M could not be ±1; when S > and 
A/ = 0, all pairs of fl^t, S, -1; 1, 0) , | 3 , S, 1; 0, 1)) form 
dark states with c.' = cj . Therefore, this kind of 
dark state is 



|5, M = Q) d = -=(\g, S, -1; 1,0) + |«?, 5, 1; 0, 1)) (16) 

Another kind of solution is = when S = N. 
In this case, there is not atomic excited state with to- 
tal spin S — N + 1, so that only the second equa- 
tion in (fl4|) is required. Thus, when S = A, all 
pairs of (\g,N,M - 1; 1, 0) , \g,N,M + 1;0, 1>) form 
dark states with ^/(A - M)(N - M + l)c+ =Ar ' M_1 = 



Thus far, it is easy to find that a dark state consist of 
two states with the same total angular momentum along 
quantize axis, Mh, and the same total spin S. The two 
states are \g,S,M - 1;1,0) and \g,S,M + 1; 0, 1), and 
their probability amplitude satisfy the equations 



N n S,M-l 

a s o 



S+1,1,M,-1 



S,M-1 Nr<S.M+l 
C + "S^S+l,!,!*,! 1 



S.M+l 







°S °S-1,1,M,-1°+ 



S,M-1 h NriS,M+l 



S.M+l 



J s w s-i,i,jwr,i l 



(14) 



Existence of solution for equation (|T4|) is the determine 
of coefficient matrix equation to zero, which gives 



y/(N + M)(N 

state is 



M + l)c 



S=N,M+1 



This kind of dark 



^/2(M 2 + TV + A 2 ) 



t/(A-M)(A-M + 1) 
v/2(Af 2 + 7V + A 2 ) 



g,N,M+l;0,l) 



(17) 



AT/JV /f~iS,M—l f~,S,M+l _ y^S,M+l ^~,a,m-i \ 
a S°S l S+l,l,M,-l°S-l,l,M,l u S+l,l,M,l u S-l,l,jV,-li 

2M(2S + 1) 



,/4S(2S- 1)(S+ l)(2S + 3) 



= 



(15) 



In equation (fT6|) and (|17p , the time factor in include into 
the state itself, because it only depend on S. 

In the above analysis, it is assumed that magnetic 
field is zero. However, it is impossible to have abso- 
lute zero magnetic field environment. Therefore, we con- 
sider how small the magnetic field is required to achieve 
dark states. Because the dark states consist of two states 
with atomic Z-component of angular momentum differ- 
ence of 2h, the Zeeman split is 2gB, where B is magnetic 
field, g = -0.70MHz/G is Zeeman split factor [20]. The 
energy difference between state with S + 1 and S — 1 
total spin is J (AS + 2). For a system with N = 1000 
atoms and uj tra p — 1000-ffz optical trap frequency, J is 
estimated to be 1Hz approximately, if Zeeman split is 
much smaller than energy split due to total spin, .i.e. 
2gB << J (AS + 2), the Zeeman split can be neglected 
and the dark state exist; otherwise, the dark states is 
damaged. It can be seen that magnetic field is require 
to be smaller than 10 _6 G so as to mean the requirement 
of dark states. Since magnetic field on the surface of the 
earth is about 0.6G, it need very accurate experiment set 
up the counteract this magnetic field. 

In conclusion, the Hamiltonian of dressed spinor BEC 
regardless of atomic and photon decay is given. The 
Hamiltonian is diagonalized under bases of eigenstates 
of non-coupling part of Hamiltonian of spinor BEC and 
optical cavity mode. Only the matrix elements that 
conserve total angular momentum along quantize axis 
is nonzero. As a result, the dressed state with total 
angular momentum along quantize axis Mh consist of 
three kind of states in non-coupling bases: |e, S', M; 0, 0), 
\g,S,M - 1;1,0) and \g,S,M + 1; 0, 1). Dark state 
is a state that oscillate between two kind of states, 



4 



\g, S,M — 1; 1, 0) and \g, S,M + 1; 0, 1), while is not ex- Acknowledgments 

cited to the states |e, S, M; 0, 0). It is found that there 

are two species of dark states, whose wave function is 

given by equation (|16p and equation (jTTJ) , respectively. 

It is obvious from the two wave functions that the dark 

state is entangle state of spinor BEC \g, S, M ± 1) and 

optical cavity mode \fi = ±). This entanglement could 
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